Thermodynamics of electron-hole liquids in graphene 
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The impact of renormalization of the electron spectrum on the chemical potential, heat capacity, 
and oscillating magnetic moment is studied. The cases of low and high temperatures are considered. 
At low temperatures, doped graphene behaves as the usual Fermi liquids with the power tempera- 
ture laws for thermodynamic properties. However, at high temperatures and relatively low carrier 
concentrations, it exhibits the collective electron-holes features: the chemical potential tends to its 
value in the undoped case going with the temperature to the charge neutrality point. Simultane- 
ously, the electron contribution into the heat capacity tends to the constant value, as in the case of 
the Boltzmann statistics. 
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I. INTRODUCTION 

Optic and magneto-optic experiments with graphene 
layers have been successfully interpreted^ so far in a 
scheme of massless relativistic particles with a conical 
energy spectrum 
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£ s (p) = =RP 



(1) 



where v is the constant velocity parameter in two bands, 
s = 1, 2, near the K and K' points in the Brillouin zone. 
In pure graphene, the chemical potential is situated at 
the charge neutrality point e = 0. However, it can have a 
nonzero value because of doping or under a gate voltage. 
Thus, the chemical potential is determined by the total 
number of carriers (difference of electrons in the upper 
band and holes in the low band) 



N = AS / |/( £ - M )-/( £ + M )| 



d 2 p 
(2irhf 



(2) 



where / (e — /i) is the Fermi function, S is the surface 
of the graphene layer, and the factor 4 takes the valley 
and spin degeneracy into account. The integration is 
performed over e > 0, the chemical potential is positive 
for electrons and negative for holes. At the fixed N, this 
condition determines the dependence p-(T), shown in Fig. 
[TJfor a relatively low electron concentration. 

For the conical spectrum, Eq. ([1]), the ratio between 
the kinetic and Coulomb energies has a constant value 
independently of the carrier concentration and the prob- 
lem of the phase electron-dielectric transition becomes 
undefined. It was recently discovered 2 in studying of 
the Shubnikov-de Haas oscillations that electron-electron 
interactions are very important for low carrier concen- 
trations, p — > 0. While the electron concentration de- 
creases from 10 12 to 10 9 cm~ 2 , the velocity parameter v 
grows by three times from its ordinary value 1.05 xlO 8 
cm/s. The logariphmic renormalization of the velocity 
for the linear electron dispersion was found by Abrikosov 
and Beneslavsky^ in the three-dimensional case and in 
Refs4ri for two-dimensional graphene. Notice, that no 
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FIG. 1: Chemical potential versus temperature for the carrier 
concentration 10 10 cm~ 2 ; the exact solution to Eq. @ is 
shown by the solid line, the asymptotes are for low, Eq. Q, 
and high, Eq. ©, temperatures (dashed and dashed-dotted 
lines, correspondingly); the renormalization is not included. 



phase transition was revealed even at the lowest carrier 
concentration. We can conclude that Coulomb interac- 
tions do not create any gap in the spectrum. 

The rcnormalized electron dispersion can be written in 
the form 



e s (p) = Tvp[l +g\n(p /p)] 



(3) 



where g = e 2 /8irhve is the dimensionless electron- 
electron interaction and e ~ 2.5 describes an effect of a 
substrate and self-screening in graphene, po ~ 0.5 x 10 8 
cm -1 is the cutoff parameter—. In Fig. [5J we consider a 



screening effect on the chemical potential at the carrier which gives the chemical potential 
concentration n = 10 12 cm~ 2 . 



Equation (j3|) is written in the linear approximation in 
g\n(po/p) < 1. Because the logarithm is assumed to be 
large, the condition g <C 1 has to be fulfilled, and we 
suppose this condition in what followed. In this article, 
we consider the impact of the renormalization on ther- 
modynamic properties of graphene such as the chemical 
potential, heat capacity, and magnetic moment. 



II. TEMPERATURE DEPENDENCE OF THE 
CHEMICAL POTENTIAL 



■k no(hv) 
41n2 f~ 



■[l + 2g\n(p v/2T)} 



(6) 



We see the inverse temperature dependence of the chemi- 
cal potential, as a collective effect in electron-hole liquids. 
The renormalization term, correcting the temperature 
dependence, is presented here explicitly and illustrated 
in Fig. H 



For low (^ 3> T) and high (/i <C T) temperatures, the 
analytical expressions for /i(T) can be obtained from Eq. 
(J2|) with the renormalization taken into account. 

For low temperatures, it is convenient to differentiate 
Eq. ([2]) with respect the temperature, using 



dfje-fijT)) 
dT 



£ — \l 

T 



dfi 
dT 



df(e-n) 



ik 



Here, we have a sharp function of (e — fi). Therefore, in 
the integrand, the momentum 

p = e[l - ghi(pov/e)]/v 

should be expand near e = fj. in powers of (e — (i) , which 
gives a factor proportional to T after the integration. For 
instance, we get in the case of electron doping 



= 



df(e-n) 



dc 



dpi (e — nY 



x[l 



dT T 

2gln(p v/^)]de, 



where we do not differentiate the logarithm because of 
the condition g <C 1. Integrating, one finds 



dfx 

dT 



3ef 



-T. 



(4) 



where we denote ef = £*(T = 0), positive for electrons 
and negative for holes. Let us notice that this is the 
known temperature dependence of the chemical potential 
in the degenerate Fermi system at low temperatures. We 
emphasize that the Fermi energy ef is determined indeed 
by the carrier concentration 



P F 1 ( £ f 
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[l-2 ff ln( P o«/|e F |)] ) (5) 



© by 



which introduces the renormalization in Eq. 
means of ef- 

For high temperatures, we can expand the integrand 
in Eq. ([2]) in /i. Introducing the new variable x = e/2T, 
we get the integral 



N 



4j4gT f 00 l-2g\n(p v/2Tx) 
n(hv) 2 J cosh 2 x 



xdx . 



III. HEAT CAPACITY 

Now we consider the electron contribution in the heat 
capacity. The energy of carriers 



E = XSl £ |/( £ - M )-/( £ + M )|J!P 

'0 (27Tft) 



(7) 



differs from the carrier concentrations, Eq. ^, only by 
the additional factor e in the integrand. Therefore, we 
can follow the same procedure. 

For low temperatures, T <C £f, the carrier heat capac- 
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FIG. 2: Chemical potential versus temperature for the car- 
rier concentration 10 12 cm -2 ; the renormalization is included 
(solid line), the chempotential for noninteracting electrons is 
shown in the dashed lines; the cutoff parameter po = 0.5 x 10 8 
cm -1 , the dielectric constant e = 2.5. 



ity in the case of electron doping writes as 
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1 — 2g\n.(pov/fj,)]d£ 
[l-2gln(p v/n)]. 



Using Eq. Q, we have 



4 e) = 2 -^-T[l-2gH Po v/\s F \)] 
3 (nv) 

in both cases of the electron or hole doping. 

For high temperatures, T 3> [m, one can perform the 
expansion of the energy, Eq. (|7|) , in the first order of fi 



E = 



it {hv)' 

2ttS\p\ 



df(e) 
de 



[1 — 2g\n(p v/e)]de 



3(hv) 

Using Eq. (J6]), we find 
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Thus, we see that the renormalization modifies the heat 
capacity at low temperatures, i.e., in the degenerate 
statistics. At high temperatures, the heat capacity pos- 
sesses the constant value and does not reveal any renor- 
malization at least to a first approximation in gln(po/p)- 



IV. MAGNETIC SUSCEPTIBILITY 

Magnetic susceptibility is determined by the depen- 
dence of the thermodynamic potential on the magnetic 
field 



n(B) 



2eBTS 
irhc 



EK 1 



in terms of the electron dispersion e sn for two bands s = 
1, 2 with the Landau number n = 0, 1, 2.. We neglect the 
spin splitting of the levels in comparison with the large 
Landau splitting in graphene. 

Oscillations of the magnetic moment in the semi- 
classical region can be found applying the Poison formula 
to the thermodynamic potential 



k=£0 J ° 



+ e t 



where the contributions of two bands are written explic- 
itly. The integraton by parts gives 
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For the semi-classical region, we use the quantization rule 
in the Bohr-Zommerfeld form 



2nn 



cA{e) 
efiB 



with the aria enclosed by the electron trajectory for the 
energy e in the momentum space 

A(e) = 7T (^) 2 (1 - 2g\n(p v/e) 

according to Eq. ([3]). 

The main contribution in the integral Q comes from 
the vicinity of the point e = ±/i(T = 0) = ±ef for 
the positive and negative sp, correspondingly. Expand- 
ing the exponent in the integrand near that points and 
integrating, one finds 

2eBTS ^ 1 sm[kcA(e F )/ehB] 

*" '~ nhc 2-^ksmh(2ir 2 kc\m(s F )\T/ehB)' 

where m(e) = ■£ — -p^ is the cyclotron mass. In calculat- 
ing of the magnetic moment we can derivative only the 
rapid factor in the argument of sin with respect B: 



u lUl .*Z*Lx 
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cos[kcA(e f) / 'efiB] 
smh[2n 2 kc\m(e F )\T/ehB} 



(10) 



where the carrier concentration no = A(e p) / (irh) 2 . This 
is the standard Lifshiz-Kosevich formula used in Ref^ for 
the interpretation of experimental data concerning the 
velocity renormalization. There are two important fea- 
tures: first, the aria A(e F ) and the effective mass m(e F ) 
should be taken at the renormalized Fermi energy cor- 
responding to the carrier concentration and, second, the 
factor in front of the sum differs from the 3d case since 
the integration over p z is absent now. 

It is interesting to compare the amplitude of oscilla- 
tions with the monotonic part of the magnetic moment, 
Ref£A 
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Thus, we see that the ration of the oscillating and mono- 
tonic parts of the magnetic moment has the order 



\M/M \ ~ 127rn 



cosh 2 (e jF /2T) 



cT 



evBJ sinh[27r 2 c|TO(£ F )|T/e?LB] 



To observe the oscillations, the argument of sinh has to be 
small or at least on the order of unity. Then, the mono- 
tonic part of the magnetic moment can be observable only 
at relatively high temperatures, \Mq/M\ ~ exp(— e F /T). 



V. CONCLUSIONS 

It should be emphasize that such a transport prop- 
erty as the electronic conductivity is not sensitive to the 
electron-electron interaction since the conductivity does 
not depend indeed on the velocity parameter v. The 
renormalization of the electron spectrum due to Coulomb 
interactions in graphene is noticeable in thermodynamic 
properties especially at low temperatures and for small 
carrier concentrations hq < 10 10 cm -2 , as can be seen 
from Eqs. ((6]), (|HJl, and (fT0|) . However, the interest- 
ing temperature dependences /i ~ T^ 1 for the chemical 
potential at high temperatures T 3> fj,, appears indepen- 
dently of electron-electron interactions. The detection of 



the renormalization requires the high accuracy in exper- 
iments because the renormalization can be concealed by 
increasing of the velocity parameter v [see Eq. ([5J and 
Fig. [5]. 
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